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SUMMARY 


Equations are derived for the stress distributions caused by three 
types of loading on infinitely long, circular, semimonocoque cylinders 
with flexible rings. These external loads can be said to lie in the 
sheet-stringer surface of the cylinder; that is, they have no radial 
components. The results, given as formulas for the stringer loads and 
shear flows in the shell, are of use in the stress analysis of circular 
semimonocoque cylinders with cutouts. In order to facilitate this appli 
cation, the formulas can be used to construct tables of influence coef- 
ficients giving stringer loads and shear flows in the neighborhood of 
each particular external load due to a unit magnitude of that load. 


INTRODUCTION 


An airplane of semimonocoque construction usually has various 
openings for doors, cargo hatches, landing-gear mechanisms, and other 
purposes. The structural design of the airplane near such openings 
requires knowledge of the stress distribution about cutouts in semi- 
monocoque shells. Stress analysis of semimonocoque structures with 
cutouts is beyond the scope of engineering beam theory because stress 
concentrations occur in the neighborhood of the cutouts. It is for 
this reason that more powerful methods of analysis are sought. 

A method for the stress analysis of sheet-stringer panels with cut- 
outs has been presented in reference 1. The method is based on the 
idea introduced by Cicala (ref. 2) that the effect of the cutout can be 
reproduced by superposing certain perturbation stress states on the 
stresses which would occur in the structure without a cutout. Three 
types of unit perturbation solutions were obtained, and it was shown 
how these solutions could be used to handle cutouts in sheet- stringer 
panels under axial loads. 



2 


WACA TN 3199 


It is evident that this approach can be extended to the analysis of 
stresses around cutouts in circular semimonocoque cylinders. In order 
to make the extension, however, perturbation solutions for this type of 
structure must be developed. Accordingly, the purpose of this report is 
to derive the perturbation solutions for an infinitely long circular 
cylinder of semimonocoque construction. These solutions are the stress 
distributions due to three types of loads analogous to those discussed 
in reference 1. The particular applied perturbation loads for this 
structure are shown in figure 1. In the analysis for the stresses due 
to these loads, bending flexibility of the rings in their own planes is 
taken into consideration. 

The use of these perturbation solutions in the stress analysis of 
circular semimonocoque cylinders with cutouts is described in detail in 
reference 3* 


SYMBOIB 


A cross-sectional area of stringers 

A n - 3B5^ - 1 + cos nS 


= Aiif n (i) 
2L sin Jjf- 



B _ E _A_ r£ 

G bt l 2 

B n = 3B5^ + 2(l - cos nS) 


arc distance between stringers, R5 


. _ *ll£n(i) 


J ni 


2L sin ^ 


(n > 2) 


ar6 

IL^b 
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En = 


2(3n - 1) 

r „ 2 


y rn 


A n 


ij 


(rm + n)^ [(rm + n)^ - l] 

coefficient in trigonometric series for 6 qj 
Y oung's modulus of elasticity 

tangential force on ring i uniformly distributed between 
stringer j and stringer j + 1 

coefficient in trigonometric series for stringer loads 

shear modulus of elasticity 


*“n(i) 

G 

H^(n,0) = D-rn cos(rm + n)0 

I * =—00 


H2(n,0) = > (-l) r D rn sin(rm + n)0 


I 

i 


k, l, r,s 
L 

M(i,0) 

m 

n 


!•=— 00 

moment of inertia of ring cross section 

longitudinal index, indicating stations (where rings are 
located) and bays 

has the value 1 when r is an integer and has the value 0 
when r is not an integer 

circumferential index, indicating stringers and panel rows 
integers 

distance between rings 
bending moment in ring i 

total number of stringers in the shell, m ^ 3 
index of the terms in a trigonometric series 
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p 

external concentrated force in the longitudinal direction 
applied to a stringer at its intersection with a ring 

p ij 

stringer load in stringer j at ring station i 

Q 

external shearing force per unit length applied about a 
shear panel 

q ij 

shear flow in shear panel (i,j) 

R 

radius to middle surface of sheet 

S 

external force in the longitudinal direction uniformly 
distributed along that portion of a stringer which lies 
between adjacent rings 

00 

D rn 2 

r=-< 

X) 

T(i,0) 

thrust in ring i 

t 

thickness of sheet 

j U ' 

total stress energy- 

V(i,0) 

transverse shear in ring i 

a ln> a 2n^l 

1 

a 3n ,a kn J 

r arbitrary constants 

) ^ m2 
4+3 

sin2 2* 

Pn = 3 + 

d 

12CS n 

± _ 3 B5 2 

2 sin2 Sp- 

7n = -2 + 

12CS n 

Aii 

second central difference in the i- direction or longitudinal 
direction, that is, A i;L g(i) = g(i + l) - 2g(i) + g(i - l) 
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rs 


5n * 


A ln' A 2n 


central angle between stringers , 2n/m 

Kronecker delta, takes the value 1 when r = s and takes 
the value 0 when r f s 


quantities defined immediately following equation (21) 
angular coordinate for rings 


X n = | cos" 1 


ifen + 1 N 

f ,3 

A 3 , 

; - 7n 


(Dn ^ l) 


x cosh"-*- 
2 




(Dn < 1) 


— ^ OAcVl - 1 

2 


tn = o cosh ' 


Pn - 1 + Mn + 1' 


" ? n 


BASIC ASSUMPTIONS 


A typical portion of a structure of the type discussed in this 
report is illustrated in figure 2. It consists of a thin-walled circu- 
lar cylinder stiffened by stringers in the longitudinal direction and 
by rings in the circumferential direction. The structure is infinite 
in length. 

The analysis in this report is based on the following assumptions 
regarding the properties of the structure: 

(a) The stringers are uniform and equally spaced around the shell, 
and the sheet is of constant thickness. 

(b) The stringers carry only direct stress. The sheet takes only 
shear stress which is constant within each bay; thus, the stringer 
stresses vary linearly between rings. 

(c) The rings have a finite bending stiffness in their own planes, 
but they do not restrain longitudinal displacements of the stringers. 
The ring bending is inextensional. 
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(d) The difference between the radius of the neutral axis of the 
rings and the radius of the middle surface of the sheet is negligible. 


ANALYTICAL APPROACH 


Equations for the stress distributions arising from the three per- 
turbation loads illustrated in figure 1 are derived in this report. The 
perturbation load shown in figure l(a) is a concentrated force in the 
longitudinal direction applied to one stringer of the shell at a ring 
station. In figure l(b) the perturbation load shown is a force in the 
longitudinal direction uniformly distributed along the portion of one 
stringer of the shell between adjacent rings. Finally, the perturba- 
tion load shown in figure l(c) consists of forces uniformly distributed 
on the portions of stringers and rings which border one shear panel of 
the shell, the forces acting in such a manner as to cause pure shear 
in that panel. These three loads will be called, respectively, the 
concentrated, the distributed, and the shear perturbation load. 

A detailed explanation of the use of these perturbation loads in 
the stress analysis of circular semimonocoque cylinders with cutouts is 
presented in reference 3« Very briefly, the idea is as follows: The 

cylinder is analyzed by any applicable method with the assumption that 
the cutout does not exist. Then a system of perturbation loads is 
applied to the shell without the cutout, and these loads are adjusted 
in magnitude to force the boundary conditions at the cutout to be 
satisfied. Stresses are obtained by superposing the perturbation stress 
distributions on the stresses which would exist in the structure without 
a cutout. The method is sufficiently versatile to handle the effects of 
reinforcement about the cutout. 

The perturbation solutions will be obtained by use of the principle 
of minimum complementary energy. This principle states that, among all 
possible stress distributions in the structure which satisfy equilibrium 
and the boundary conditions on stress, the distribution that most nearly 
satisfies compatibility is the one which minimizes the complementary . 
energy it* where 


n* = Internal energy - Work done by Surface stresses acting through the 

prescribed displacements 


( 1 ) 


Since displacements are not prescribed anywhere on the structure, the 
second term on the right-hand side of equation (l) is omitted. The com- 
plementary energy becomes the internal energy or stress energy of the 
structure. 
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In writing the equation for the stress energy, the following factors 
are considered: the energy of axial distortion of the stringers, the 

shear energy in the sheet, and the bending energy of the rings in their 
own planes. In figure 2 are shown the index systems used to designate 
stringers, rings, bays, and panel rows. Note that the intersection of 
ring i and stringer j occurs at the lower left-hand corner of shear 
panel (i,j). Each of the perturbation loads is shown in its positive 
sense in figure 1. Stringer loads are taken as positive in tension. 

Shear flows are positive as shown in figure 3. Ring bending moments, 
shear, and thrusts are placed on the ring element in figure 3 in the 
positive sense. With this index notation, the stress energy in the 
structure can be expressed as 



where the integration over the length of a stringer between adjacent rings 
has been carried out. 

In the analysis to follow, stringer loads are expressed in the form 
of a finite trigonometric series. Then, by using the equations of statics, 
the shear flows and ring bending moments are written in terms of the 
coefficients of this trigonometric series. The expression for stress 
energy, equation (2), is minimized with respect to the coefficients of 
the trigonometric series for stringer loads j then, the expressions for 
the stringer loads , shear flows , and ring bending moments are substituted 
into the resulting equation. This process yields a fourth-order finite- 
difference equation which can be solved for these trigonometric coeffi- 
cients. The solution is then substituted back into the original expres- 
sions for stringer loads, shear flows, and ring moments to yield the 
desired distributions . 

For convenience in application, the significant equations are col- 
lected in appendix A. 
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PERTURBATION LOAD SOLUTIONS 
Concentrated Perturbation Load 


Expression for stringer loads .- The concentrated perturbation load 
is shown in figure 1(a); let P represent the magnitude of this load. 
Since the structure is uniform and infinitely long, half of the load 
goes into the portion of the structure to the right of the ring station 
where the load is applied .(ring i = 0), and half goes to the left of 
this ring station. Therefore, it can be seen from figure l(a) that, 
because of symmetry. 


Pij - "P-i,,j 

r— 1 

All 

(5a) 


(i i 0) 

(3b) 

l(i,0) = -M(-:i,0) 

(i So) 

(3c) 


Consider the right half of the structure, including the ring at 
i = 0. The concentrated perturbation load gives rise to stringer loads 
which are circumferentially symmetric about stringer j = 0 (see fig. 1(a)). 
Thus the stringer load distribution can be' represented by a series of the 
form ' 


S or 
2 


m-1 

2 


■ ij 


= V f n (i)cos nj& 


n=0 


GO 


where the notation 
from n = 0 to n 


m Qr m-1 
2 2 

^ means that the summation is carried over n 

n=cT 

= S if m is even and to n = m - - if m is odd. 


Evaluation of fo(i);. 
is multiplied hy cos 2j5 
procedure yields 

Pi a cos l = 
j=0 1J 



fj_(i)^ and f n (0).- Suppose that equation (if) 
and summed over j from. 0 to m - 1. This 


S§ or m - 1 


m-l 


f n(i) ^ cos nj6 cos Ij 5 
n=0 j=0 
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The sum over j on the right-hand side is, for 0 % n < S 0 S l S 


m-1 


cos nj6 cos IjB = 0 


«j=0 


(l i n) 


+ &n0 + 5 n ,|) (l “ n) 


Thus the coefficients of the trigonometric series in equation (4) are 

m-1 


fn (l) »- S_ 

mf 1 + &n0 + 5 n ^mj j=0 


Pij cos nj5 


(5) 


It is desirable first of all to evaluate those f n (i)'s which can 
be found from consideration of the boundary conditions and of the over- 
all equilibrium of the cylinder. Consider the equations of statics for 
the cylinder as a whole. Satisfaction of equilibrium in the longitudinal 
direction requires that the sum of the stringer loads at any ring sta- 
tion i must equal one-half of the applied load P. This condition is 
expressed as 




For n = 0, equation ( 5 ) gives 


V 1 ) 



(6) 


Moment equilibrium gives two equations, one of which is automatically 
satisfied because of the symmetry of the stringer load distribution around 
the cylinder. The other moment equation is 

m-1 



10 


NACA TN 3199 


For n = 1, equation ( 5 ) is 


m-1 


f l(i) =I^IPij cos =| 

J ^ 


(7) 


On substituting the values of fg(i) and f]_(i) given in equations ( 6 ) 
and ( 7 ), respectively, into equation (4), there results 


- a + I 008 J s + 


m m -1 
2 or ~ 


n =2 


f n (i)cos nj5 


(8) 


Consider now the boundary condition at ring station i = 0. The 
stringer loads here are 


p 0 j _ 2 S 0 j 


and substitution of this expression into equation ( 5 ) yields 


f n(0) 


P 



+ &n0 + 8 

n 



(0 < n < |) (9) 


The equations of equilibrium and the boundary condition at i = 0 
have been used to obtain certain of the coefficients of the trigonometric 
series for stringer loads • The remainder of the coefficients are found 
by use of the principle of minimum complementary energy, and this is the 
next step in' the solution. 


Expressions for shear flows and— ri-ng- bending -moments In- order- to- 
use the principle of minimum complementary energy, the shear flows and 
ring bending moments must be found in terms of the trigonometric coeffi- 
cients fn(i). Shear flows are determined by the consideration of the 
equations of statics of a portion of any stringer j between two adjacent 
rings i and i + 1. The forces on this free body are shown in the sketch 
below: 


q. .L 
ij 


P U 


^j-l 1 


p i+l>j 


L 
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Equilibrium of these forces requires that 

Pi+1,0 " Pij + (Wj " ^i,j-i) L = 0 

Substitution of equation (8) into equation (10) yields 

n- 
2 


m or m-1 


" q i,j-l 


n=2 


( 10 ) 


|Pn( i+1 ) " f n (iT]cos nj6 (ll) 


In order to find q^j, replace j with a dummy index k and sum both 
sides of this equation over k from k=l to k = that is, write 


( q ik " q i,k-l) " L 


| orjS|l 


ff n (i+l) - f n (ij~| ]>I cos nk& 
n=2 k=l 


When the indicated summations over k have been carried out, the fol- 
lowing equation is obtained: 


S or !ti 

2 _ 2 


9ij - 1i0 


n=2 


[f n (i+l) - f n (i)J 


sin n^j + ^5 1 


2 sin ^ 


The term q^Q can be found from the condition that the total torque on 
the section is zero. The resulting expression for shear flows is 


| or-^- fn(i+l) _ fn(i) 


m-1 




n=2 


2L sin 3^- 


sin 


n (j - |): 


(12) 


Bending moments are caused in each ring by a tangential loading 
which develops because of the difference in shear flow in the sheet on 
either side of the ring. The tangential load on ring i has the value 

m-1 


2^2 




' = " S 


n=2 


ALi^nU) 

2L sin 2^- 


sm 


n (j * 


(13) 
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In appendix B, this load is applied to a circular ring and the bending 
moment in the ring is derived. This procedure results in the following 
moment in ring i (see eq. (B9)): 


M(i,0) = " 2 YL 2 A ii f n(i)Hl(n,0) (l4) 

n=2 

where 


Hi(n,0) = 


cos(rm + n)0 


r “°° (rm + n)^ [(rm + n)^ - l] 


The sign convention for the moment is illustrated in figure 3; the con- 
vention for measuring angle 0 is shown in figure 4(a). 

Energy analysis .- Now the stringer loads , shear flows, and ring 
bending moments have been expressed in terms of the coefficients f n (i). 
The stringer loads are given in equation (8), the shear flows in equa- 
tion (12), and the bending moments in equation (l4). . These' equations 
are used in the minimization of the stress energy of the cylinder with 
respect to fn(i). 

By virtue of the symmetry properties of this problem expressed in 
equations (3a), (3b), and (3c), the energy in the structure to the left 
of ring i = 0 is the same as the energy to the right of ring i = 0. 
Thus, equation (2) becomes 



Note that M(O,0) is identically zero because there is no difference in 
shear flow across ring i = 0 and, therefore, no tangential load acts on 
this ring. 
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Minimization of the stress energy with respect to fn(i) results 
in the following equation: 


dU 


d f n (i) 


= 0 


m-1 

t. 

j=0 


L 

6ae 


(pi + i , j + **Pij + + lr(%j + 


/ 


dq. . 


a n-i,j' 


■*2>t 


'0 


R_ 

El 




M(i-1,0) 


dM(i-l,0) 

df n (i) 


10 


(15) 


The coefficients fg(i) and fp(i) are known already for all values of i, 
and f n (o) is known for 0 -g n < S. Thus equation (15) needs only to he 

considered for i > 1 and n 2? 2. The expressions for the stringer loads, 
shear flows, and ring bending moments are substituted into equation ( 15 ). 
Then the following definite sums are needed (these can be obtained by 
the procedure outlined in ref. 4): 


m -1 


cos njS = 0 


0=0 


(0 < n < m) ( 16 ) 


and for the integers n and 2 restricted to the range H n<- and 


1 < 2 < — 
1 = - 2 # 


m -1 


cos ZjS cos nj& = 0 


0=0 


C + Vf ) 


and 


^i 

J 


> sin 2 (j + ^jb sin n^j + = 0 


- if 1 + 8 


(2 0 n) 
(2 = n) 

V 

(2 0 n) 


m ] ( ^ ~ n ) 

’ 2 . 


> (IT) 


H 18) 
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The following definite integral, which is derived in appendix C, is also 
needed: 


where 



H 1 (Z, 0 )H 1 (n, 0 )d 0 



(19) 




(rm + n) 1 * - [(rm + n) 2 - 1_ 


and where n and l are restricted to 2 ^ n ^ | and 2 < 2 ^ E. a 
closed form of S n is presented in appendix C but the series form con- 
verges so rapidly that it is usually more convenient than the closed 
form for use in calculations. 


After substitution of the expressions for stringer loads, shear 
flows, and ring moments into equation (15), the use of these definite 
sums (l6), (17), and (l8), and definite integral (19) results in the 
following equations which express the condition of minimum stress 
energy: 

For i = 1 , 

f n ( 3 ) + 2r n f n (2) + (2p n - l)f n (l) + 2( 7n + l)f n (0) = 0 (20a) 

and, for i ^ 2, 

f n (i+2) + 2 7n f n (i+l) + 2p n f n (i) + 2 7n fn(i-l) t f n (i-2) = 0 (20b) 

t 

where 


1 _ 3 B6 g 
2 sin2 VL 
2 


?n = " 2 + 


12 CSn 
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4 + 3 


Pn = 5 + 


B6^ 

sin 2 ^ 


12CS 


n 


B _ E _A_ 

G bt l 2 


C = 


AR° 

IL^b 


Solution of finite-difference equation ,- Equation (20b) is a fourth- 
order finite-difference equation with constant coefficients. (Note that 
the symbol i represents the index of the rings and bays and should not 
be confused with the usual notation for i£l which never appears in 
this report.) Equation (20b) corresponds exactly with equation (24) of 
reference 5 . The general solution is presented on pages 23 to 2 6 of 
reference 5 and on pages 28 and 29 of reference 6. It may be written as 


fn(i) = (±e * n ) j^lrAln(i) + a 2n A 2n( i j] 

(±e^Y ja 3 n A ln ( i ) + 04^(1 j] 


(n ^ 2) (21) 


where the upper sign is used when y n < 0 and the lower sign when 
7 n > 0. The A's are as follows: 


2(Pn ~ 1) 

For Dn = — 2 — > 1, 

7n 


A ln(i) = cos iV 


A 2n( i ) = sin i*n 
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where 


' X 



cos 


-1 



1 





for D n < 1, 


A Ln(i) = cosh iX n 


where 


A^Ci) = sinh iX n 




i cosh"-*- 



2 




and, for 1^ = 1, 


Am = 1 

A 2n = i 

In the inverse trigonometric and hyperbolic functions, the principal 
values are used. The argument \|r n of the exponential function is given 
by the positive branch of 


^n 


- cosh - "'' 


3n - 1 



At a large longitudinal distance from the applied load, the stringer 
loads should approach the elementary distribution given by the first two 
terms of equation (8); consequently, for n ^ 2, f n (i) — >0 as i— ><». 
The first term on the right-hand side of equation (2l) satisfies this 
condition; however, the second term does not and, hence, must be omitted. 
The solutions , then, that are compatible with the boundary conditions at 
infinity are: 7 


f n(i) - Cn f^liAlnC 1 ) + a 2n A 2n^ i ^] 


(n > 2) (22) 
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where 

tn - ±e-*n 

Now the arbitrary constants and. are determined. The first 
a]_Q is obtained immediately. Substitution of i = 0 into equation (22) 
and use of equation (9) to evaluate f n (o) yields 


(n ^ 2) (23) 


Substitution of equations (23) and (22) into the boundary equation (20a) 
yields 


fn(o) - a iji " 


m 


-C 1 


m 


= - 


®ln + * 1) 


9 2n “C 1 + 8 n,§) 


where 


©sn = Zn A snO) + 2 7n £n 2 Asn(2) + (2p n - l)^ 8 n (l) (s = 1,2) 


The solution for the concentrated perturbation load is now complete 
since the coefficients f n (i) are completely defined and may be substituted 
into equation (8) to give the stringer loads. The shear flows can be 
found from equation (12); however, once the stringer loads are known, 
it is simpler to calculate the shear flows by the use of the equations 
of statics. Because of symmetry, the shear flows in shear panels adjacent 
to stringer j = 0 are given by 


. _ wo - P 1+1,0 

q i0 ^ 1,-1 2L 


All the other shear flows can be found by the use of equation (10). If 
desired, the moment distribution in the rings can be computed from 
equation (ik) and the thrust and transverse shear in the rings can be 
found from the formulas given in appendix B. 
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Distributed Perturbation Load. 

Expression for stringer loads .- The distributed perturbation load 
is shown in figure 1(b); let S represent the magnitude of the total 
force distributed along stringer j = 0 between rings i = 0 and i = 1. 
From figure l(b) it is seen that 


p ij = -p-i+1, j ' 

(i S 1) 

( 24 a) 


(i i 1) 

( 24 b) 

M(i, 0 ) .= -M(-i+l, 0 ) 

(1 S 1) 

( 24 c) 

At ring station i = 1 and to the right of this 

ring, the 

stringer 


loads can be represented by a trigonometric series of exactly the same 
form as equation (4) 


8 - 




m-l 

2 


f n (i)cos nj& 


n=0 


( 25 ) 


except now i = 1, and the coefficients f n (i) are different from those 
obtained for the preceding case of the concentrated load. 


Evaluation of fo(i) and f]_.(i).- As in the preceding case, the 
first two coefficients fo(i) and f^(i) can be obtained from the equa- 
tions of statics, and the results are the same as before. Equation (25) 
becomes 


p i0 = S + I COS J6 + 


m ^ 

2 or 


m-l 

2 


f n (i)cos nj& (i = l) (26) 


n=2 


With the concentrated perturbation load, all the coefficients f n ( 0 ) 
were easily found because the stringer load distribution at ring sta- 
tion i = 0 was known. Here no such distribution is known. In order to 
determine the boundary condition at bay i = 0, the effect of the dis- 
tributed perturbation load on the equilibrium of portions of stringers 
in this bay must be investigated. 

Expressions for shear flows and ring bending moments .- Away from 
bay i = 0 the shear flows and ring bending moments are of the same form 
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as for the concentrated load. The following expression for the shear 
flows is obtained by use of equation (10): 


m m-1 
2 ° r ~P~ 




f n (i+l) - f n (i) 


n=2 


2L sin 


n6 

T 


sin n^j + ^jb (i > l) (27) 


The ring bending moments are obtained in appendix B as 


where 


M(i,0) 


B or ti 
2 . 2 


n=2 


2rtL 


Aiif n (i)Hi(n,0) 


(i ^ 2) (28) 


OO 

Hi(n,0) = ]T 


cos(rm + n)0 


]?=— OO 


(rm + n ) 2 (rm + n ) 2 - lj 


The applied force in bay i = 0 may be written as S6oj* Consider, 
now, the equilibrium of a portion of any stringer j between ring i = 0 
and ring i = 1. The forces on this free body are shown in the sketch 
below: 


qoj L / 

:w / 




p 0j 






90,4-1 l 


■^ p lj 


Equilibrium of these forces requires that 

p lj " p 0j + (<K>J “ 9 o,j-l) L " S 5 0 j = 0 
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Because of the antisymmetry property expressed in equation (2^-a), the 
equilibrium equation becomes 

^lj + (^0j " ^J-lX ' ^Oj = 0 (29) 


It is convenient, now, to expand the Kronecker delta 8 q^ in a 
finite trigonometric series, 


& 0j = 


m 

2 or 


m-l 

2 


d n cos njS 


n=0 


(30) 


Multiplying through by cos ZJ8 and summing over j from 0 to m-l, 
yields the trigonometric coefficients d n . The result is 


^n 


2 


m ^- + &n0 + 


(31) 


Substitution of the expression for stringer loads, equation (26), and 
the trigonometric expansion for 6ojj equation (30), into the equilibrium 
equation (29) yields 


^0j " q O,j-l 


2 2L_— 


[sd n - 2f n (l)j cos njS 


n=2 


In order to find q^j, this equation can be treated in the same manner 
as equation (ll); that is, replace j by a dummy index k, sum from 
k = 1 to k = j, and then use the condition that the total torque on 
a cross section in bay i = 0 must be zero. This procedure results in 
the following expression for the shear flows in bay i = 0: 



n=2 


| Sd n - f n (l) 

L sin 3* 

' 2 


sin 


* §> 


(32) 


The expression for the bending moment in rings i = 1 and i = 0 is 
yet to be found, as this expression differs from that for the moment in 
the rest of the rings given in equation (28). The moment in ring i = 0 
is the same in magnitude as that in. ring i = 1 but opposite in sign. 
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The tangential loading on ring i = 1 is given by 


S or Shi 

2 _ — 2 f n ( 2 ) - 3fn(l) + Sd n 

q u “ = " Z_ " sin 

n =2 


2L sin ^ 


n ( 3 * ±> 


By analogy with equations (13) and (l4), then, the bending moment in 
ring i = 1 can be written as 


M(l,0) = - 


“ or Shi 
2 _ 2 


R2 m 


n =2 2 jtL 


[£n(2) - 3f n (D + Sd^Hifn,#) (33) 


All the stringer loads, shear flows, and ring bending moments have 
now been expressed in terms of the coefficients f n (i). The stringer 
loads are given in equation ( 26 ), the shear flows in equations ( 27 ) and 
(32), and the ring moments in equations ( 28 ) and (33). It re ma ins to 
substitute these expressions into the equation obtained from minimization 
of the stress energy of the cylinder with respect to f n (i). 

Energy analysis .- By virtue of the symmetry properties in this 
problem, given in equations (2k), the energy in the structure to the 
right of bay i = 0 equals the energy to the left of this bay. Equa- 
tion ( 2 ) for the stress energy can be written 






RSL 2 
2Gt qi <3_ 



R 

2EI 


M^(i, 0 )d 0 


(3k) 


Minimization of the stress energy with respect to f n (i) results in the 
following equations: 
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dU 

dfn(D 



+ 




6 f n (l) 


+ Sl/a*. i2i i. + 

2GtV lj df n (l) 


^OJ 




JR 

El 


N(2,0) 


dM(2,0) 

3f n a) 


+ 


M(1,0) 


c>M(l,0) 


d 0 


and 


dfn(D 



i+ 1,0 


+ 4p 


10 


+ Pi-l, 





(35) 




_R 

El 


M(i+ 1 , 0 ) 


dM(i+l,0) 

df n (i) 


+ 


M(i,0) 


aM(i,0) 

*n<*> 


+ 


M(i-1,0) 


dM(i-l, 0 ) 


d 0 


(i ^ 2) (36) 


Note that equation ( 36 ) is the same as equation ( 15 ), except here it is 
valid only for i ^ 2 . 

The stringer loads , shear flows, and ring moments are substituted 
into equations ( 35 ) and ( 36 ), and then the definite suras and definite 
integral derived in the preceding section are used to simplify these 
equations • After simplification, the following equations result: 

For i = 1, 


/p n - 47 n - 2\ 

f n (3) + (2r n " l)f n (2) + 2 (Pn " 7 n )f n (l) = Sd n^— T~ ) (37a) 
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for i = 2, 

fn(M + 27 n f n (3) + 2p n f n (2) + (2 7n - .l)f n (l) = -Sd n (3Tb) 

and, for i = 3> 

f n (i+2) + 27 n f n (i+l) + 2p n f n (i) + 27 n f n (i-l) + f n (i-2) = 0 (37c) 

Solution of finite-difference equation .- Equation (37 c ) is the same 
as equation ( 20b ) ; therefore, the solution to equation (37 c ) is 

f n(i) = Cn 1 lo-^lnC 1 ) + a 2n A 2n( i )] ( n = 2 ) (38) 

which is the same as equation (22) except for the values of the arbitrary 
constants cv\_n and These constants are found by the substitution of 

the solution (38) into equations (37 a ) and (37b). This procedure yields 
two simultaneous algebraic equations in aq n and an ^ their solu- 

tion gives 


a ln = 


Pn “ ^7 n - 2 

r 2n + n 2n 

3 


2S 


n ln r 2n “ r ln a 2n 


n ( 1 + V| 


a 2n = - 


Pn " ^7 n - 2 


fl ln + F ln 3 

n ln r 2n " r ln^2n 



where d n , the coefficient in the trigonometric series for the Kronecker 
delta 5oj- > ^ as been replaced by its value as given in equation (3l)> 
and where the ft's and f's are given by 


&sn _ Cn^sn(3) + (27 n ~ l)£ n ^A3n(2) + 2(3n ~ 7n)Cn^sn(^) ( s _ 1>2) 

(39a) 

Isn = Cn^sn(^) + 27 n £ n ^Asn(3) + 2p n ^ n ^A sn (2) + 

(2 7n - l)^ sn (l) (s = 1,2) (39b) 
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The coefficients f n (i) are now defined for the distributed perturbation 
load and may be substituted into equation (2 6) to give the stringer loads. 
The shear flows can be found from equations (27) and ( 32 ), but, again, 
once the stringer loads are known shear flows can easily be found by 
use of the equations of statics. The shear flow in the panels adjacent 
to stringer j = 0 can be found by considering symmetry: 


In bay i = 0 


^0 " ~%,-l ~ 


S ~ ^10 
2 L 


and, outside of bay i = 0, 


< liO = 


-^-l = 


P i0 “ P i+1,0 


2 L 


(i >*1) 


The other shear flows are found from equation (lO), as before. If 
desired, the ring moments can be obtained from equations (28) and (33) 
and the ring thrust and transverse shear can be found from the equa- 
tions given in appendix B. 


Shear Perturbation Load 

Expression for stringer loads .- The shear perturbation load is shown 
in figure 1 (c). The magnitude of the load per unit length applied along 
the stringers and rings bordering shear panel (0,0) will be represented 
by Q. From figure l(c) it is seen that the longitudinal symmetry prop- 
erties in this case are the same as those for the case of the distributed 
perturbation load given by equations ( 24 ). 

The shear perturbation load is self-equilibrating and gives rise 
to stringer loads which are antisymmetric about panel row j = 0 . For 
i > 1, the stringer loads may be represented by 


Pij = f n (i)sin n(j - ( 40 ) 

where the coefficients f n (i) are different from those in the two pre- 
ceding cases. The term corresponding to n = 1 vanishes because it 
represents an elementary bending stringer-load distribution, and the 
shear perturbation load does not require this distribution for overall 
equilibrium. 
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Expressions for shear flows' and ring lending moments .- None' of the 
coefficients f n (i) in the trigonometric series (kO) can be found from 
the equations of statics. Fur thermore , the boundary condition at 
bay i* = 0 must be determined from a consideration of the effect that the 
shear perturbation load has' on the equilibrium of the. portions of stringers 
in bay i = 0 and on the bending moment in the rings bounding this bay. 

Thus the energy approach must be used immediately and the first step in 
this approach is to write the shear flows and ring moments in terms of 
f n (i), the coefficients of the trigonometric series for the stringer 
loads, equation (^0). 

Outside of bay i = 0, the satisfaction of the equations of statics 
for the portions of stringers between adjacent rings yields equation (10), 
the same as in the two preceding cases. Substituting equation (4o) for 
the stringer loads into the equilibrium equation (10) and following the 
same procedure used to obtain equation (12) yields the expression for 
the shear flows due to the shear perturbation load: 




S or 

2 2 f n ( i+ D . f n(i) . 

> cos nj5 

n=2 2L sin — 


(i ^ 1) 0+1) 


The tangential loadings on the rings to the right of ring i = 1 are 




S or 51=i 

2 sr — 2 Aiifn(i) 


n=2 2L sin 


nS 


cos nj& 


In appendix B this load is applied to a circular ring and the following 
expression for the moment in the ring is obtained (see eq. (B13)): 


m or m-1 


M(i,0) = - 2 X 2 Sr Aiif n (i)H 2 (n,0) 
n=2 2ltL 


(i ^ 2) (42) 


where 


E 


w 

: 2 (n,0) = (-i) r 


sin(rm + n)0 


r=-oo 


(rm + n)^ |(rm + n)^ - l] 


The convention for measuring angle 0 here is a little different than 
before and is illustrated in figure 4(b). 
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Now, the shear flows in bay i = 0 and the bending moments in the 
rings bordering bay i =' 0 must be found. Consider the shear, flows in 
this central bay. The shear perturbation loading applied at bay i = 0 
may be written C&oj* Then tlie forces on the portion of any stringer j 
between ring i. = 0 and ring i = 1 are as shown in the sketch below: 


/ 


^0j L 


qOjL 


/ 




V 


/ 


•Oj 




r : 

tr — ^ 

q 0,j-i L 

\ 

\ 




L — * 



p lj 


Ojj-l 1 


Equilibrium of these forces requires that 

p lj " p 0j + (^Oj - ( 30,j-]) L + Q ( 5 0,j-1 - 6 0j) L = 0 


Because of the antisymmetry property, equation (24a), the equation of 
equilibrium becomes l 

^lj + (^j " + < 5 0 ,j-l " & 0 j) L = 0 (^3) 

The substitution of the stringer loads 7 equation (40), into 'the equilibrium 
equation (43), and the introduction of the trigonometric expansion for the 
Kronecker delta Sqj, equation (30), yields the following equation: 


m rsv. m ~l 
2 


_q O,j-l = "I ZL 2 f n (l)sin n^j - - Qd L [cos(j - l )8 - cos j&] 


S or . St! 
2 . 2 


n =2 


Qd n [cos n(j - 1)5 - cos njb] 
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Now q^j can be found by replacing j with a dummy index k, summing 

over k from k=l to k = j , and using the condition that the torque 
on a cross section within bay i = 0 balances the applied torque. This 
procedure results in the following equation for the shear flow. in the 
central bay: ' . ...... 


IS or — 
2 _ _ 2 


qoj = Q^o + 0^1 cos + 


^n(l) 


iL sin 


nS 


+ Q£ n ]cos njS 


(MO 


Consider the bending moment in rings i = 1 and i = 0.. The moment 
in ring i =0 is identical in magnitude to the moment in ring i = 1 but 
of opposite sign. The tangential loading per unit arc length on the 
portion of ring i = 1 between stringer j and stringer j + 1 is illus- 
trated in the sketch: 


\ 


r 

> 



1 












£ stringer j + 1 


£ stringer j 


Adding these tangential loads and introducing the series expansions for 
q^j and q^j and for 6 qj yields the total load per unit arc length 
on ring i = 1: 


^lj - qoj + ^0j " 



n=2 


fn(2) - 3f n (l) 
2L sin ^ 


cos nj& 


By analogy with equations (13) and (l4) the bending moment in ring i = 1 
is 


M(l,0) 


m 

2 



|3|F n (2) - 3f n (l[|H2(n,0) 


(^5) 
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Expressions for stringer loads, shear flows, and ring moments have 
been written in terms of the coefficients f n (i). The stringer loads are 
given in equation (40), the shear flows in equations (4l) and (44), and 
ring moments in equations (42) and (45). These expressions are ready 
to be substituted into the equation which results from minimizing the 
stress energy with respect to fn(i). 

Energy analysis . - Because the longitudinal symmetry relations which 
exist for the distributed perturbation load, equations (24), also exist in 
the case of the shear perturbation load, the stress-energy expression used 
in the distributed-load problem can be used here. The expressions obtained 
on minimizing this stress energy, equations (35) and (36), are also appli- 
cable here. So the stringer loads, shear flows, and ring moments just 
derived are substituted into equations (35) and (36). At this stage in 
the two preceding cases, certain definite sums and a definite integral 
were introduced to simplify the equations. A similar procedure is 
followed here. ' 

The definite sums which are of interest are 

m-1 / \ 

y~ sin nfj - i)5 = 0 
j=0 ^ 


and for the integers n and Z restricted to the range 1 <L n ^ 2 an a 

ls.sj ■ ' " 2 


and 


m-1 



cos Zj6 cos nj6 = 


0 



m-1 

sin Z 
j=0 



sin n 



(Z / n) 
(Z = n) 


(Z / n) 


_ m 
2 


( 1 + 


n 



( l = n) 
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is 

(Z / n) 

= «S n ^L + B n ^mj ( Z = n) 

where n and Z are restricted to 2 ^ n ^ S and 2 Z < jS. . After 
simplification the following equations result: 

For i = 1, 

fn(3) 


+ (27 n - l)f n (2) + 2(p n - 7 n )f n (i) = “2L#L n 


'Pn- - ^7n - ll\ 


sin 


nS 


The definite integral, which is derived in appendix C, 
- 2 * 


/ H 2 (Z,0)H 2 (n,0)d0 = 0 


for- i = 2, 


(46a) 


f nOO + 27 n f n (3) +.2p n f n (2) + (27 n - l)f n (l) = 0 (46b) 

and, for i ^ 3> 


fn(i + 2) + 27 n f n (i+l) + 2p n f n (i) + 27 n f n (i-l) + f n (i-2) = 0 (46c) 


Solution of finite-difference equation !- Equation (46c) is the same 
finite-difference equation for which the solution is written in the two 
preceding sections. Substitution of this solution, equation (38), into 
equations (46a) and (46b) gives two simultaneous algebraic equations 
for aqn and a 2n , the arbitrary constants. Solution of this system 
yields 


2n 


P n - 4 7n - 11 


. ’ nS 
sin — 
2 


4QL 


fi ln r 2n " r ln n 2n 


m ( 1 + 6 n,§) 


a ln = - 
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0-211 = 


n Pn - ^7n - 11 nS 

r ln 3 sin — 

fi ln r 2 n " r ln^ 2 n 


AQL 


m (i + 6 n; |) 


The ft's and r's in this case are precisely the same as in the pre- 
ceding case of the distributed perturbation loadj ft sn is given by equa- 
tion (39a) and ^sn by equation (39h). 

With the coefficients f n (i) known for the shear perturbation load/ 
the stringer loads are obtained from equation (AO) and the shear flows 
can be found from equations (1+1 ) and (Ah) . For panel row j = 0, the shear 
flow equations become 


and 


<liO 


5 — ^ fn(i+l) - f n (i) 


n=2 


2L sin Sp. 


(i ^ 1) 


3Q , 


m m-1 

2 or T 


qoo = — 


m 


n=2 


f n(D 


2Q 


| L sin f- m(l + 


When the shear flows in panel row j = 0 are. known, it is simpler' to 
compute the remainder of the shear flows by use of the equations;, of 
statics rather than equations (hi) and (Ah). In shear panels (0,l) 
and (0,-l) adjacent to the loaded panel, the shear flow is given by 


2p lr + QL 

%1 “ < 30,-1 " 9oo £ 

All the other shear flows are found by use of equation (10). If desired, 
the ring bending moments can be found from equations (A2) and (A5) and 
the ring thrust and transverse shear can be calculated from the formulas 
given in appendix B. 
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LIMITING CASE OF RIGID RINGS 


If the ring bending stiffness is allowed to increase indefinitely, 
the rings approach complete rigidity in bending, the parameter C approaches 
zero, and a considerable simplification results. For this limiting case, 
equations ( 20 ) for the concentrated perturbation load reduce to 

f n (i+i) - 2 2 tt f n (i) + f n (i-i) =0 (i £ 1 ) ( 47 ) 

A n 

where 


A n = 3B6^ - 1 + cos nS 


B n = 3B6^ + 2(l - cos nS) 


This can be shown easily by multiplying equations (20) through by C 
and allowing C to approach 0. Equation (47) is a second-order finite- 
difference equation with constant coefficients. The same equation, 
together with its general solution, is given in reference 6 , page 31 . 

The solution compatible with the boundary conditions at infinity can be 
written as 


fn(i) 


o«(±e-' n 


) l 


(48) 


where 


cosh A n = 


Bn 

An 


and where the upper sign is taken when 


A n 


< 0. 


A n > 0 


and the lower sign when 


The arbitrary constant <x n is determined by evaluating the solution, 
equation (48), for i = 0 and introducing the value of f n ( 0 ) given in 
equation ( 9 ). The result is identical to given in equation (23) 


P 


On = 


C + 6 1 ,§) 


(n > 2 ) 
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Equations (8) and (12), the expressions for stringer loads and shear 
flows, respectively, used before in the case of the concentrated pertur- 
bation load are still valid. The substitution into these expressions of 
the solution (48) with the constant 0 ^ as found above yields the stringer 
loads and shear flows due to a concentrated perturbation load when the 
rings can be considered rigid. 

For the case of the distributed perturbation load, equations (37) 
reduce in the limit to 

■ (-A n )f n (2) + (2B n + An)fn(l) = 3B5 2 Sd n 

fn(i+l) - 2 f n (i) + fn(i-D =0 (i > 2) 

The arbitrary constant an in the solution (48) is 


6B5 2 S 

A»( ± e-*n t 7 ) "( r + ^l) 

For the shear perturbation load, equations (46) reduce to 


(-A n )f n (2) + (2B n + A n )f n (l) = -6LQd n B6 2 sin ^ 

f n (i+l) - 2 f n (i) + f n (i-l) = 0 
A n 

The solution is again equation (48) and becomes 

12B5 2 sin f- ul 

On = - — — 

A n (*e‘ >ta + l) m ( 1 + 6 n,| ) 


(i S 2) 
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CONCLUDING REMARKS 


Solutions have been obtained for the stress distributions caused 
by three types of loading on infinitely long, circular, semimonocoque 
cylinders with flexible rings. The solutions are adaptable to the con- 
struction of tables of influence coefficients which give stringer loads 
and shear flows in the neighborhood of the perturbation load due to a 
unit magnitude of that load. The perturbation loads are the tools of a 
method described in NACA Technical Note 3200 for the stress analysis of 
circular semimonocoque cylinders with cutouts. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 22, 195^* 
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APPENDIX A 

SUMMARY OF SIGNIFICANT EQJATIONS 


The formulas and parameters required for computing the stress dis- 
tribution due to concentrated, distributed, and shear perturbation loads 
are collected in this appendix for reference. 


Stringer Loads 

Concentrated perturbation load (see fig. 1(a)): 

S or 2-1 

p 2 __ 2 

p ij = ail + m COS + f n(i)cos nj& 

n=2 

where P is the applied load. 

Distributed perturbation load (see fig. 1(b)): 

m m-1 

2°i_- 


a a _ _ 

Pii = — + - COS + J> f n (i)cos nj& 

0 2m m n=2 


(i ^ 0) 


(i £ 1) 


v/^here S is the total applied load. 

Shear perturbation load (see fig. 1(c)): 

m-1 

— or 

, _ 2 ^ 
ij 


- or 

2 2 

f n (i)sin 


n=2 


(i ^ 1) 


Shear Flows 

Concentrated perturbation load (see fig. 1(a)): 
For the shear panels in panel row j = 0, 

_ p i0 " P i+1,0 

5 10 - — — 
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and, for the remainder of the shear panels. 


_ _ p ij " p i+l;j . _ 

l 


(j ^ 1) 


Distributed perturbation load (see fig. 1(b)): 
For the shear panel (0,0), 


^0 


S - 


^10 


2L 


for the remainder of the shear panels in panel row j = 0, 


_ Pio - Pi+1,0 

qi 0 . — il — 


(i i 1) 


and, for all other shear panels, 


nd - Pij ~ r Pi * i,j + 


(j ^ 1) 


Shear perturbation load (see fig. l(c)): 
For the panel about which the load is applied, 


= 5Q 

m 


m ^ m-l,_ 

5 or 


<300 = tr + 


n=2 


fn(l) 


2Q 


L sin ^ mi 


( P + # “«i) 


for the remainder of the shear panels in row j = 0, 


9i0 


21 nr Sill 

2^-2 f n (i+l) - f n (i) 


2L sin 


n=2 


(i £ 1) 
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for the shear panel (0,l), 


^01 " ^00 


2P n + 03 L 


for the remainder of the shear panels in panel row j = 1, 


Pil " Pi+1,1 
q il £ + q i0 


and, for all other shear panels. 


q ij = ^ uL - ~ - £ i+1,j + q i,j-l 


Evaluation of the Trigonometric Coefficients fn(D 


Basic parameters: 


for Flexible Rings 


-d E _A_ r£ 

B " G bt l 2 


Auxiliary parameters: 


4 + 3 


3 B5*^ 

^ sin^ — 
2 
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Discriminating parameter: 


Dn = 


2(p n - 1) 



Trigonometric coefficients: 

^n(^) = ^n |^ln^ln(^) + c-PryA pr/ ijJ 

where 


Cn 



Aln(i) = cos iX n 
= 1 

= cosh ix n 


A2n(i) = sin iX n 
= i 

= sinh iX n 



(n ^ 2) 


(Dn > 1) 
(D n = 1) 

( D n < 1) 
(Dn > 1) 

(Dn = 1) 
(D n < 1) 

(Dn > 1) 

(Dn < 1) 



38 


NACA TN 3199 


Arbitrary constants for concentrated perturbation load: 


a ln = 


"f- + 6 n,f ) 


c^n = - 


®ln '*’ 2 (7 n + l) P 


®2n 


m(i + »„«) 


where P is the applied load and 


0, 


*sn £n ^sn(3) + 2 7 n £n ^sn( 2 ) + ( 2 3 n - ljl^Ag^l) (s = 1,2) 


Arbitrary constants for distributed perturbation load: 


a ln = 


P _ 3n " ^7n - 2 . n 

l 2n y— + ^2n 

fl ln^2n ” r ln^2n 


2S 


+ 6 n,|) 


B ln * r ta 


P n - ^7n - 2 


°2n = 


2S 


fl ln r 2n " r ln a 2n 


m 


c + 


where S is the total applied load and 

ft sn = £ n 3A sn(3) + ( 2 7 n " lK n 2A sn(2) + 2(p n - 7 n K n A sn(l) 


(s = 1,2) 


f’sn ^n + 2 7n£n 3 ^sn(3) + 2 3n£n Asn( 2 ) + ( 2 7n _ n£ n Ag n (l) 


(s = 1,2) 
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Arbitrary constants for shear perturbation load: 


P n - 4 7n - 11 


sin 


am = - 


n ln r 2n ‘ r ln n 2n 


f + 5 n,|) 


Pn - **?n - 11 


am = 


ft ln r 2n " r ln n 2n m 


(1 + 6 m\ 

\ n ’2J 


where Q is the applied load per unit length. 


Evaluation of the Trigonometric Coefficients f n (i) 


Basic parameter: 


for Rigid Rings 


E A_ R^ 
G bt L 2 


Auxiliary parameters: 


A n = 3B£> - 1 + cos n6 


B n = 3B6 2 + 2(l - cos nS) 


An = cosh' 


-1 Pn 


Trigonometric coefficients: 


f n (i) = Onf— e" An 
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Arbitrary constant for concentrated perturbation load: 


^ = 


P 

"f - + 6 n,|) 


Arbitrary constant for distributed perturbation load: 


6B6 2 S 



Arbitrary constant for shear perturbation load: 
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APPENDIX B 

BENDING MOMENT, AXIAL THRUST, AND TRANSVERSE SHEAR IN RINGS 


Expressions will be developed for the bending moment, axial thrust, 
and transverse shear in a circular ring under tangential loads such as 
those which arise from the differences in shear flow across a ring in a 
circular semimonocoque cylinder. 

Two cases must be considered: One case occurs with the concentrated 

and distributed perturbation loads, where the ring loading is anti- 
symmetric about stringer j = 0. The other case occurs with the shear 
perturbation load, where the ring loading is symmetric about panel 
row j = 0. 


Concentrated and distributed perturbation loads .- For the concentrated 
and distributed perturbation loads, the tangential loading on ring i has 
been written in the form of a finite trigonometric series (see eq. (13)) 






(Bl) 


where 


J nl 


ALi^n(i) 
2L sin 


(n £ 2) 


This ring load has a stepwise variation around the ring, being constant 
between stringers and having jump discontinuities at the stringers. The 
limitation that n ^ 2 ensures that the ring is in equilibrium. 


The procedure will be to expand each term of the series (Bl) in an 
infinite Fourier series in the variable 0. For each harmonic of the 
Fourier series, that is, for a continuous sinusoidal tangential force 
distribution on the ring, the moment, thrust, and shear in the ring are 
easily found (see ref. 5, p. 33, for example). On the basis of inexten- 
sional deformation and the neglect of transverse shear distortions, the 
results are as follows: If the tangential load on ring i is given by 


cos n0 + b n ^. sin n0 


(n £ 2) 
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then the moment, thrust, and shear in this ring sire, respectively, 


. ,, R 2 , _ 

*^(i,0) = “ a ni - 7 -z 77 sin n0 + *> n i 


R £ 


n(n2 - 1) ‘ ni n (n 2 - l) 


cos n0 


“* H — R 

T n(i>0) = -a ni — 7 n sin n0 + b ni — n cos n0 )> (B2) 


n 2 - 1 


n 2 - 1 


£ p 

V n (i,0) = a^i — cos n0 + b ni — sin n0 


n 2 - 1 


Figure 3 shows the sign convention used in writing equations (B2). 

Consider, now, one term of the series (Bl). To expand this term 
in a Fourier series, write 


-m 


sin 


00 

n (j + §)& = 2 _ ( c r)ni sin r 0 
\ / r=n 


(B3) 


where the ( c r) n i' s ane the Fourier coefficients. It is obvious that 

the first harmonic which will occur in the Fourier series in equation (B3) 
must be that for which r = n. The other harmonics, then, will be added 
to this to build up the step shape of the loading function. The con- 
vention for measuring angle 0 in this - case is illustrated in figure 4(a). 
The index j can be thought of as a function of 0, that is: when 
0 < 0 < 5 , j = 0; when 6 < 0 < 2&, j = 1; and so forth. 

In order to carry out the expansion, equation (B3) is multiplied 
through by sin Z0 and integrated from 0 to 2* 

ro-il. n(j+l)S v r>2jt °° 

2 _ / 'bni sin n/j + ±J& sin 20 d0 = / (c r ) ni sin r0 sin 20 d0 

j=0 J j6 ' ‘ J 0 r=n 

After integration, the right-hand side of this equation becomes 
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by virtue of the orthogonality of the trigonometric functions. The left- 
hand side becomes 


2 sin nfcl 

1 b ni 2_ sin + 2 P Sin Z ( j + 2 r 


on carrying out the integration. From reference' 4 it can be shown that 


n (j + sin i(j + |)5 = | (-1) m J z _ n - (-1) m J z+n 


Thus the Fourier coefficients are given by 


( C Z ) • 
v 1 ni 


sin 


- 1 ) J 2-n - 


/ \ HI 

(- 1 ) Jl 


The nth term of the tangential loading on the ring is 


b n j_ sin 


Z+n . Z6 

s in -tt 


n (j + i)& = S b ni 21 (-1) m J z - (-1) m Jz+n — sin Z0 

v ‘ * i=n -zr ~j r l 


By use of the properties of J r this summation can be rewritten 


b ni sin n (j + |) 5 = f b ni 21 C” 1 ) 


sin(rm 4- n)% 


sin(rm + n)0 


(- 1 )’ 

r=l 


c 

sin(rm - n)^ 


rm - n 


sin(rm - n)0 
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On expansion by- the sum and difference formulas of trigonometry and with 
the use of the fact that m& = it is found that 


sin(rm + n)| = (-l) r sin ^ 


> 


(B 6) 


sin(rm - n)^ = (-l) r+1 sin 

When equations (b 6) are substituted into equation (B 5 ), the following 
relationship results: 


b ni sin n^j + = S b ni sin ^ 


Z 


sin(rm + n)0 ^ sin(rm - n)0 


_q rm + n r= i rm - n 


sin(rm + n)0 


m . . n& V slI H rm + : 

= z ^ni sin > : 

n d r==oo rm + n 


(BT) 


From the first of equations (B 2 ) it is seen that if the tangential 
loading on the ring is given by the right-hand side of equation (B7) 
then the bending moment in that ring is 


Mn(i, 0 ) = R 2 5 b ni sin ^ H^n^) (b8) 

where 


00 

Hi(n, 0 ) = ]T 


cos(rm + n)0 


r =-oo 


(rm + n ) 2 |(rm + n ) 2 - lj 


Equation (b 8) gives the bending moment in a ring which carries a 
tangential load distributed according to one term of the series of 
equation (Bl). When the ring is loaded by the sum of such stepwise 
terms, as in equation (Bl), then the moment is given by a sum of terms 
like (B8). So the bending moment in ring i is 


M(i,0) 


m 

2 


2 


m-1 
- 2 


n=2 


R 2 - b 


ni 


sin Hq(n, 0 ) 


(B9) 
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For completeness, the expressions for axial thrust and transverse shear 
can be written in a similar manner - 


T(i,0) = 


5 or 

2 r — 2 


n=2 


R ® *ni sin f Ki(n,0) 


V(i>0) = 


•2 or 

2 .r — 2 


n=2 


R S b ni sin 2& L x (n,0) 


where 


K x (n,0) 



cos(rm + n)0 

— 

(rm + n) - 1 


Ll (n,0) = Z 8 JL n .<™ + n)0 


r= '“ (rm + n) [Tnn + n) 2 - l] 


Shear perturbation load .- In the case of the shear perturbation 
load, the tangential loading on ring i is given by the finite trigo- 
nometric series 


where 


F ij " *lj - 9i_i,j - 



n=2 


(BIO) 


Aijfn(i) 

2L sin 


(n = 2) 


Equation (BIO) can be treated in a manner analogous to the handling of 
equation (Bl). That is, each term of the series in equation (BIO) can 
be expanded in a Fourier series. Then the moment, thrust, and shear in 
the ring are written immediately. 
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Analogous to equation (B3 ) , write 


a n i cos nj5 = >_ (c r ) ni cos r 0 


r=n 


where, now, the angle (f> is as shown in figure 4(b). If both 
tion (Bll) are multiplied by cos Z0 and integrated from 0 
results for the Fourier coefficients: 


(c 1) 


ni 


2 sin ^ - 
2 _ 

Til 


m-1 

cos nj5 cos Ijt 

j=0 


It can be shown (see ref. 4) that 

m-1 
x — 

J 


cos nj6 cos IJ8 - § + Jjj) 
J \ m m / 


so the nth term of the tangential loading on the ring is 


^ - N sin 

ani cos nj6 = - ani 2_( J Z-n + J Z+n) cos l f 


l=n \~ It) 1 


This summation becomes 


a nl cos nJS = f ^ sin § J (-l) r 


r=-« 


rm + n 


which corresponds to equation (B7). Then the bending moment 

ffl or 

M(i,?l) = - 21 2 = %I sin f. H 2 (n,jZl) 

n=2 rt d 


Similarly, thrust and shear are 


| or 2ti 
2 T 2 


T(i,0) = - 21 R S ani sin if K 2 (n,0) 


n=2 


(Bll) 


sides of equa 
to 2jt, there 


(B12) 


is 


(B13) 
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V(i,0) = 


m m-1 

2 °L_— 


® ~ a ni si n ^ ^2 ( n j 0) 


, ,, V - , .r sin(rm + n)0 

H 2 (n, 0 ) = (- 1 ) - — 7 = — =r 

r=-« (rm + n)^ Km + n)^ - 1 


K 2 (n,0) = 


sin(rm + n)0 
(rm + n)^ - 1 


Ie(n,0) = £ (-l) r ^l™.+ n)0 - 

r 00 (rm + n) |(rm + n)^ - lj 
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APPENDIX C 


EVALUATION OF DEFINITE INTEC2ZAIS 


In order to minimize the stress energy it is necessary to investigate 
the following definite integrals: 

p2n r2.it ° p . °° _ 

/ Hi(n,0)H 1 (Z,0)d0 = / D ra cos(rm + n)0 \ D s i cos(sm + l)0 d0 

0 u 0 r=-°° s=-°o 


(Cl) 


and 


P 2n r 2n r 2 - r 

/ H2(n,0)H2(l,0)d0 = / ^ (-1) D rn sin (rm + n)0 \ (-l) s D s j sin(sm + z)0 d0 

J 0 Jo r=-oo s=-oo 


(C2). 


where 


®rn - 


(rm + n)^ [(rm + n)^ - l] 
and where integers n and Z are limited to the following ranges: 


2 ^ n ^ | 

2< l<- 
~ “2 


Consider the relation (Cl). The right-hand side can be written 


00 00 


h Z Z V”: 


I*=-oo S=-°° 


sm + Z)0 d0 


si 


P 2rt 

J cos(rm + n - 


r 2it 

sm - Z)0 d0 + / cos(rm + n + 
J 0 


00 00 


" 2 J D rn D s ^5 rm+ . n ^ sm+ .2 + Srm+n,-sm- l)2jt 


00 g;=— 00 
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l*-9 


Now, by virtue of the limited range of the integers n and l, the following 
relations can be written: 


^rm+n,srrB-Z ^ c 


Z-n _ ^sr^ Zn 




1 
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it is found that when n = — 

2 



Hlff, 0)d0 = 2 


Z D rffl 2 >" 2S B” 

f =_oo >2. 2 


To summarize, then, 



Hl{n,0)H;L(l,0jd0 = 0 


(Z / n) 


= S n it (1 + 5 


i m ) 

n ’2J 


(Z = n) 


Consider the relation (C2). It is handled in a manner analogous' to 
the treatment of (Cl). Equation (C2) can be written 


p2rt 


00 00 


/ H 2 (n,0)H 2 (z,0)d0 = X I(-l)\n(-l) S Ds4( 8 »n,s W r^ n; . sm . l ) 

(J x’— -00 s=-oo N 7 / 


2* 


00 00 


S?— ^ _1 ^ r (“ 1 ^ SD rn D sz(^ & sr & nZ“ 5 s,-r-l 5 _ n£ ? mV 

r=-°o s=-o° \ L * 2 .) 


For 2^n<5, 



H 2 (n,0)H 2 (Z,0)d0 


0 


(Z ^ n) 



S n « 


(Z = n) 
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If n 


m 

2 ’ 


C H 2 (§,0> 2 (f, ! «>s = £ 


m 


- (-i) r (-i)' r 'V 


r=-oo 


m 


p 


m 



- 2S m it 

2 


Thus the definite integral (C2) gives precisely the same result as (Cl) 


H 2 (n,0)H 2 (l,0)d0 =0 (l / n) 


= S n rt^l + & n ^m^ (l = n) 


The sum 



can he expressed in closed form with the aid of formula 6 . 495 , number 2, 
reference 7 * The result is 

g _ 5 ^ 2 + cos nS + 6 2 cos nS cos S - 1 5 5 sin S 

(l - cos nb)^ T - cos n& 4 ( cos _ cos g)^ 4 cos ng _ cos g 



However, the series form of Sn, because of its rapid convergence, may 
be more convenient than the closed form for use in computation. 
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Figure I Perturbation loads . 



Figure 2.- Portion of typical cylinder. 
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Figure 3.- Positive sense of quantities used in analysis. 




(a) For concentrated and (b) For shear perturbation load . 

distributed perturbation loads 


Figure 4.- Conventions for angular coordinate <f> . 
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